Effect of disorder on the vortex-lattice melting transition 
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We use a three dimensional stacked triangular network of Josephson junctions as a model for 
the study of vortex structure in the mixed state of high T c superconductors. We show that the 
addition of disorder destroys the first order melting transition occurring for clean samples. The 
melting transition splits in two different (continuous) transitions, ocurring at temperatures Ti and 
T p (> Ti). At Ti the perpendicular-to-field superconductivity is lost, and at T p the parallel-to- field 
superconductivity is lost. These results agree well with recent experiments in Y BaCuO. 



Although the existence of a melting transition of the 
vortex structure of high T c superconductors is now ac- 
cepted, the nature of the transition and its dependence 
on different parameters is not completely clear in spite 
of the big deal of theoretical and experimental work per- 
formed 0]. Experimentally, it is known that some ex- 
ternal parameters can have strong effect on the behavior 
of the system. Some of these parameters are the mag- 
netic field, the anisotropy, and the quenched disorder. It 
was found that in some cases (low fields and weak disor- 
der) the transition is first order || whereas in other cases 
(high fields or strong disorder) it is continuous J3|,|j . 

Theoretical studies on the nature of the transition fo- 
cus on the problem of whether it is a single transition at 
which superconducting coherence is lost along all direc- 
tions of the sample, or two successive transitions where 
perpendicular and parallel to field coherence is lost in 
two sequential steps [61. In addition, different possibili- 
ties for the transition such as first order, second order, or 
crossovers have been proposed. Theoretical treatments 
have been based on models that can usually account only 
for part of the experimental data. The results depend on 
the parameters used and sometimes they look contradic- 
tory §. 

In this paper we fill part of the gaps between theory 
and experiment by studying the three dimensional (3D) 
Josephson junction array (JJA) model in the presence 
of an external magnetic field and with different kinds of 
disorder. We show that when temperature increases this 
model gives different results depending on the disorder: 
for low disorder we obtain a single first order transition 
where superconducting coherence is lost in all directions, 
whereas for high disorder superconducting coherence is 
lost in two successive continuous steps, first perpendicu- 
larly and then parallel to the applied magnetic field. This 
behavior is closely related to that found in experiments 

The equations for the 3D JJA model are la] 
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where if 1 (t) is the phase of the superconducting order pa- 
rameter, I c and Rq are the critical current and normal 
resistance of each junction, and 4>q is the flux quanta. 
Eq. ([j]) gives the current j" between nearest neighbors 
nodes i and i'. Here A 1 " 1 is the vector potential of the ex- 
ternal magnetic field, and i?" (t) is an uncorrelated gaus- 
sian noise which incorporates the effect of temperature. 
Eq. (||) ensures the current conservation on each node, 
and jl xt is the external current applied at node i. We 
study this model on a stack of triangular two-dimensional 
planes coupled by vertical links. 

In all the simulations shown below the external cur- 
rent j ext (which is used to calculate resistivities) is about 
~ 10~ 2 of the Josephson junction critical current in the 
corresponding direction. If we want to minimize the pin- 
ning effect of the subjacent JJA lattice we have to use a 
magnetic field H as low as possible. However, a low value 
of H increases the sample size needed to get good statis- 
tics. We have used a magnetic field H equal to 1/6 flux 
quanta per plaquette. This value of field generates an 
Abrikosov lattice which is commensurate with the subja- 
cent triangular JJA lattice and is a compromise between 
low fields and not too long computation time ||. Self- 
inductance effects are not considered -this is equivalent to 
take the magnetic penetration depth A — * 00. When we 
simulate anisotropic systems we diminish the c axis mean 
critical current of the junctions and at the same time in- 
crease the c axis elemental resistance by the same factor 

V 2 ( V 2 = (^'^ / ' {l(t) where || and _L indicate parallel 
an perpendicular to the ab plane, and (...) indicates the 
mean value on the sample). Boundary conditions (BC) 
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are used as follows: when external current is applied in 
the x direction we take open BC for the x direction and 
periodic BC for y and z directions; when external current 
is applied in the z direction we take pseudo periodic BC 
for the z direction jlO| and periodic BC for x and y di- 
rection. In all cases the temperature is measured in units 
of the mean Josephson energy of the in-plane junctions. 



resistivity for T > T m is very close to the value corre- 
sponding to decoupled planes, whereas for T < T m it is 
clearly different from that value. This indicates that an 
effective decoupling of the planes is occurring at T m . As 
we argue below, the nature of the first order transition in 
thick samples is different from that of an isolated plane. 

It has been found experimentally that disorder alters 
the previous picture. In fact, for samples with twin 
boundaries the transition in p a f, is continuous To 
study the effect of impurities in our model we allow the 
critical currents I c of each junction to vary randomly be- 
tween two fixed values I™ ln and 7™ ax . This generate pin- 
ning centers close to the junctions with low critical cur- 
rents. All the results shown in the simulations with disor- 
der correspond to D = (/™ ax - If in ) / (/™ ax + If in ) = 
0.5. 



FIG. 1. Resistivity along the ab plane for a decoupled sys- 
tem of 18 x 18 x 18 planes (open symbols), and for a three 
dimensional system with different anisotropies (full symbols, 
anisotropics rj 1 — 5, 10, 20, and 50 from right to left). 
The first order melting temperatures T m for the different 
anisotropies are indicated by arrows. Temperature is mea- 
sured in units of the Josephson energy of the in-plane junc- 
tions. 

We consider first the behavior of a system of uncoupled 
(triangular) planes when the c coupling is turned on. It 
has been found that a single plane has a weakly first or- 
der transition [jllj . Our results for the resistivity of this 
system are shown in Fig. |l| (open symbols). Within the 
numerical precision we cannot distinguish from the be- 
havior of the resistivity between a continuous transition 
and a weak first order transition occurring at T ~ 0.23 
We will refer to this type of behavior as a con- 
tinuous transition, indicating that no discontinuities arc 
observed in the p(T) curve. Both continuous (second 
order) phase transitions and crossovers are in this cat- 
egory. When the coupling is increased we clearly see a 
well defined transition temperature where the resistivity 
has a jump. This is a first order transition as found in || 
and [|l3|. In addition, the c axis resistivity (not shown) 
has also a jump at the transition temperature, indicating 
that the superconductor coherence is lost discontinuously 
and at the same temperature T m in all directions. Hys- 
teresis loops, which for the sake of clarity are not shown 
in the figure, are observed for all the cases studied but 
for the uncoupled planes. We conclude that the c axis 
coupling transforms a continuous transition into a first 
order transition 113] . Note that for high anisotropies, the 



FIG. 2. Resistivities along the ab plane (full symbols) and 
along the c direction for an isotropic 18 x 18 x 18 sample 
and different configurations of the disorder. For the clean 
sample hysteresis loops around T m upon heating and cooling 
are shown. For the disordered case the approximate values of 
Ti and T p are indicated. 

The effect of disorder depends on its spatial correla- 
tion. We consider three different possibilities: point (un- 
corrected), columnar (correlated along c axis), and pla- 
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nar (correlated in 'twin boundaries' perpendicular to the 
applied current) defects. For the uncorrelated disorder 
no correlation on the values of I c for different junctions 
are set in. For the columnar pinning all junction having 
the same ab plane coordinates have the same critical cur- 
rents, whereas there is no correlation for junctions with 
different ab plane coordinates. For the twin boundaries 
case al junctions with the same x coordinate {j ex t || %) 
have the same critical currents, whereas there is no cor- 
relation for junctions with different x coordinates. 

The results for p a b obtained for the different types of 
disorder are illustrated in Fig. || (full symbols) for an 

isotropic = (l<t)) sample. For the clean sample 

the hysteresis loop is shown, indicating again a first or- 
der transition at T = T rn . Although hysteresis in the 
resistivity is not an unequivocal indication of a first or- 
der transition, we checked that the transition is in fact 
first order by using the histogram technique Jf5|||. In 
addition other indication of loose of coherence (in par- 
ticular the helicity modulus) give results which are in 
agreement with the conclusions obtained from resistivity 
measurements. In all disordered cases the transition in 
p a b is continuous, and we denote the transition temper- 
ature as Ti. Point defects and twin boundaries decrease 
the transition temperature, whereas columnar defects in- 
crease it. This is consistent with the fact that columnar 
pinning generates the strongest pinning. 

It is instructive to study the behavior of p c with disor- 
der. In Fig. H (open symbols) we show the values of p c for 
the same disorder configurations as before. The main re- 
sult is that p c also becomes continuous in the presence of 
disorder, and the transition temperature (that we denote 
T p |^[[l0) and experimentalists T t h fl(|) is higher than the 
corresponding Ti. The temperature T p is higher than T m 
for samples with columnar defects and twin boundaries, 
whereas is lower than T m for point defects. This indi- 
cates that point defects diminish the coherence along the 
c direction, whereas the other two types of defects en- 
hance it. Note also that the absolute difference T p — Ti 
is the highest for the twinned sample. This is due to the 
fact that planar defects parallel to the vortex movement 
have opposite effects on Ti and T p : Ti diminishes because 
there are paths of easy movement for the vortices, and 
the rigidity of the ideal Abrikosov lattice has been de- 
stroyed, and T p increases because the defects are c axis 
correlated, and thus they enhance the coherence along 
the c direction. 

We found that if the value of the disorder D is smaller 
than a critical value D c (D c ~ 0.3 for isotropic samples) 
the first order transition persist and the coherence in all 
directions is lost at the melting temperature. Strong dis- 
order (D > D c ) is necessary to produce the behavior of 
Fig. | 

Our results are in complete agreement with recent re- 
sults on YBaCuO samples, where a first order transition 



is observed in clean samples and two transitions are ob- 
served with Ti < T m < T p in twinned samples R], It 
has also been shown that point defects destroy the first 
order transition H , but unfortunately pseudo DC trans- 
former or direct c axis resistivity measurements were not 
performed on these samples. 

We checked that the transitions at Tj and T p occurring 
in our triangular samples with disorder have the same 
features described in references ||JT^,|T^] for the case of 
a square lattice of junctions: the transition at Ti is a 
thermal depinning of vortex lines, which is probably re- 
lated to the large amount of disorder |Q; the transition 
at T p is a percolation phase transition of the vortex lat- 
tice perpendicularly to the applied field. Also we have 
checked that around T p , p c (T) satisfies the percolation 
scaling laws of (fj] . In addition, the structure of the lat- 
tice for Ti < T < Tp is that of a disentangled vortex 
liquid. The range T p — Ti in which this disentangled liq- 
uid exists shrinks to zero when increasing the thickness 
of the sample @Jf7[. 

When two transitions at different temperatures Ti and 
Tp occur, it is natural that both of them are continu- 
ous: at Ti the system changes from a solid to a liquid of 
c axis correlated (and thus effectively two dimensional) 
vortices. This transition is qualitatively similar to the fu- 
sion of a two dimensional system, and it is expected to be 
continuous. At T p the transition is driven by a percola- 
tion of vortex loops, which is a second order (continuous) 
phase transition. 

Experimental and theoretical evidence suggest that the 
transition is first order only when the superconducting 
coherence is lost in all directions at the same tempera- 
ture. This suggests that the first order melting in clean 
samples is a consequence of the interplay between two 
different transitions: the depinning of single vortex lines 
-which drives the crossover occurring at T^-, and the per- 
colation of vortex loops between planes -which is respon- 
sible for the appearance of dissipation along the c di- 
rection at T p . These two transitions may cooperate and 
merge onto a single one for clean samples: when perco- 
lation occurs the planes decouple, and this reduces the 
effective value of Ti. In turn, when vortices depin, vor- 
tex loops between planes are more easily generated due 
to the screening effect of mobile vortices, thus reducing 
Tp. The combination of these two effects can generate an 
instability that drives the transition first order [^9[ . 

As we said above, the merging of Ti and T p onto a sin- 
gle first order melting temperature T m is observed only 
for clean samples. However, if the thickness of the system 
is lower than a minimum value two different transitions 
at Ti and T p are observed. In particular, in our sim- 
ulations the thickness of the system has to be greater 
than ten planes on isotropic samples in order to obtain 
a first order transition. For thinner samples we observe 
two separate and continuous transitions. One can make 
a rough estimate of the dependence of the critical thick- 
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ness on field in the following way. Previous results [ pL 0|Jl 7| 
indicate that T p depends on thickness as 



ksTp ~ A/ In (L c /d), 



(3) 



with d a distance of the order of the interlayer spacing, 
L c the thickness of the sample, and A an energy which 
is the order of the energy necessary to make two near- 
est vortices touch each other. Being A proportional to 
the distance between vortices we can write A = eo<f>Q\ / ^H, 
with eo an energy scale factor related to the linear energy 
density of a vortex line. To be able to merge onto a single 
transition, TJj and T p should be of the same order. Hav- 
ing into account that Tj is rather thickness independent 
|^o| , ^7| we obtain the minimum value L m in necessary to 
have a first order transition as 



L„ 



dexp ( e 4>o/V HksT, 



(4) 



This expression shows that L m in — > oo for H — > 0, 
indicating that very thick samples are necessary to ob- 
serve the first order transition at very low fields. This 
is probably not an experimental limitation due to the 
small value of d, but it should be taken into account in 
numerical simulations at low fields. 

We have concentrated on the effect of disorder on the 
melting transition at a fixed value of the external field. 
Finding the behavior of the system when changing the 
magnetic field H is difficult because of commensurability 
effects between the vortex lattice and the mesh. However, 
we can choose to change at the same time the magnetic 
field and the discretization parameter, in such a way that 
for any real field we have 1/6 flux quanta per plaquette. 
The change in the discretization produces a change in 
the anisotropy r\ and the disorder of the system. For a 
real change in H of the form H — > H8 with 6 ~ 1, the 
anisotropy changes as r\ — > r\b x l 2 p0| . The change in 
the disorder may depend on the particular realization. 
In the simplest case of uncorrelated disorder it can be 
shown that D changes as D -> DS 1 / 2 . We conclude 
that we can study the H — T phase diagram of high- 
T c 's by analyzing the i] — D — T phase diagram along 
lines of fixed value of rj/D. Preliminary results indicate 
that the overall properties of the H — T phase diagrams 
of YBaCuO and BiSrCaCuO are well reproduced |l9| . 
In particular, from these arguments and Fig. [l] we can 
conclude that the jump in the resistivity at the first order 
melting should decrease when increasing H. 

In summary, we have presented results obtained using 
the three dimensional JJA model that reproduce very 
well the behavior of the vortex lattice in high-T c super- 
conductors. For clean samples the system has a first or- 
der melting transition at temperature T m , at which the 
superconducting coherence is lost discontinuously in all 
directions. In the presence of disorder the transition sep- 
arates in two continuous ones occurring at temperatures 



Ti and T p (> Tj). At T.- L perpendicular-to-field correla- 
tion is lost due to depinning of vortices from the pinning 
centers. At T p parallel-to-field correlation is lost due to a 
percolation transition of the vortex lattice. In addition, 
relative values of Ti, T p , and T m depend on the kind of 
disorder considered. 
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